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Abstract 

The classical Landau/Levich models of liquid- 
propellant combustion, which serve as seminal exam- 
ples of hydrodynamic instability in reactive systems, 
have been combined and extended to account for a 
dynamic dependence, absent in the original formu- 
lations, of the local burning rate on the local pres- 
sure and/or temperature fields. The resulting model 
admits an extremely rich variety of both hydrody- 
namic and reactive/diffusive instabilities that can 
be analyzed in various limiting parameter regimes. 
In the present work, a formal asymptotic analysis, 
based on the realistic smallness of the gas-to-liquid 
density ratio, is developed to investigate the com- 
bined effects of gravity, surface tension and viscos- 
ity on the hydrodynamic instability of the propa- 
gating liquid/gas interface. In particular, a com- 
posite asymptotic expression, spanning three distin- 
guished wavenumber regimes, is derived for both cel- 
lular and pulsating hydrodynamic neutral stability 
oundaries A p (k), where A p is the pressure sensi- 
tivity of the burning rate and k is the disturbance 
wavenumber. For the case of cellular (Landau) in- 
stability, the results demonstrate explicitly the sta- 
bilizing effect of gravity on long-wave disturbances, 
the stabilizing effect of viscosity and surface tension 
on short-wave perturbations, and the instability as- 
sociated with intermediate wavenumbers for critical 
negative values of A p . In the limiting case of weak 
gravity, it is shown that cellular hydrodynamic insta- 
bility in this context is a long-wave instability phe- 
nomenon, whereas at normal gravity, this instability 
is first manifested through 0(1) wavenumber distur- 
bances. It is also demonstrated that, in the large- 
wavenumber regime, surface tension and both liquid 
and gas viscosity all produce comparable stabilizing 
effects in the large-wavenumber regime, thereby pro- 
viding significant modifications to previous analyses 
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of Landau instability in which one or more of these 
effects were neglected. In contrast, the pulsating 
hydrodynamic stability boundary is found to be in- 
sensitive to gravitational and surface-tension effects, 
but is more sensitive to the effects of liquid viscosity,’ 
which is a significant stabilizing effect for 0(1) and 
higher wavenumbers. Liquid-propellant combustion 
is predicted to be stable (i.e., steady and planar) 
only for a range of negative pressure sensitivities 
that lie between the two types of hydrodynamic sta- 
bility boundaries. 


1. Introduction 

The burning of liquid propellants is a funda- 
mental combustion problem that is applicable to var- 
ious types of propulsion and energetic systems. The 
deflagration process is often rather complex, with va- 
porization and pyrolysis occurring at the liquid/gas 
interface and distributed combustion occurring ei- 
ther m the gas phase 1 or in a spray. 2 Nonetheless, 
there are realistic limiting cases in which combus- 
tion may be approximated by an overall reaction at 
the liquid/gas interface. In one such limit, the gas 
flame occurs under near-breakaway conditions, ex- 
erting little thermal or hydrodynamic influence on 
the burning propellant. In another such limit, dis- 
tributed combustion occurs in an intrusive regime, 
the reaction zone lying closer to the liquid/gas in- 
terface than the length scale of any disturbance of 
interest. Finally, the liquid propellant may simply 
undergo exothermic decomposition at the surface 
without any significant distributed combustion, such 
as appears to occur in some types of hydroxylam- 
monium nitrate (HAN)-based liquid propellants at 
low pressures. 3 Such limiting models have recently 
been formulated, 4 ' 5 thereby significantly generaliz- 
ing earlier classical models 6 - 7 that were originally 
introduced to study the hydrodynamic stability of a 
reactive liquid/gas interface. In all of these investi- 
gations, gravity appears explicitly and plays a sig- 
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nificant role, along with surface tension, viscosity, 
and, in the more recent models, certain reaction- 
rate parameters associated with the pressure and 
temperature sensitivities of the reaction itself. In 
particular, these parameters determine the stabil- 
ity of the deflagration with respect to not only clas- 
sical hydrodynamic disturbances, but also with re- 
spect to reactive/diffusive influences as well. Indeed, 
the inverse Froude number, representing the ratio of 
buoyant to inertial forces, appears explicitly in all 
of these models, and consequently, in the disper- 
sion relation that determines the neutral stability 
boundaries beyond which steady, planar burning is 
unstable to nonsteady, and/or nonplanar (cellular) 
modes of burning. 8,9 These instabilities thus lead 
to a number of interesting phenomena, such as the 
sloshing type of waves that have been observed in 
mixtures of HAN and triethanolammonium nitrate 
(TEAN) with water. 3 Although the Froude number 
was treated as an 0(1) or larger quantity in these 
studies, the limit of small inverse Froude number 
corresponding to the microgravity regime is increas- 
ingly of interest and can be treated explicitly, leading 
to various limiting forms of the models, the neutral 
stability boundaries, and, ultimately, the evolution 
equations that govern the nonlinear dynamics of the 
propagating reaction front. 

In the present work, 15,16 we formally investi- 
gate both the normal- and reduced-gravity parame- 
ter regimes in the context of a theory that exploits 
the realistic smallness of the gas-to-liquid density ra- 
tio p. The resulting analysis and results thus facili- 
tate a comparison of some of the features of hydrody- 
namic instability of liquid-propellant combustion at 
reduced gravity with the same phenomenon at nor- 
mal gravity. In addition, the small-p limit allows a 
tractable synthesis of the separate models described 
in previous analyses to simultaneously account for 
the effects of viscosity (both liquid and gas), sur- 
face tension and pressure sensitivity of the burning 
rate on this phenomenon. In contrast, the original 
analysis of Landau 6 neglected viscosity, while that 
of Levich 7 retained liquid viscosity, but neglected 
surface tension and gas viscosity. Neither of these 
classical studies accounted for any variation of the 
local burning rate due to pressure perturbations, 
which has recently been introduced into this type of 
analysis. 4,5 The latter study, in particular, presented 
a highly tractable analysis of the inviscid case, thus 
effectively generalizing Landau’s model to account 
for a pressure-dependent burning rate. By develop- 
ing a perturbation analysis in the limit of small p, 
we show how all the phenomena just described may 
be simultaneously included, for both the normal and 


reduced gravity cases. In addition, the pressure cou- 
pling in the model predicts not only the cellular type 
of hydrodynamic instability attributed to Landau, 
but also a pulsating instability associated with suf- 
ficiently negative burning-rate pressure sensitivities. 
Both forms of hydrodynamic instability are analyzed 
here, leading to explicit asymptotic representations 
of the neutral stability boundaries. 


gas/liquid interface: 



Figure 1 

Model geometry. The liquid/gas interface x 3 = 
$ 5 (xi,X2,^) separates the unbumed liquid propellant 
from the gaseous products. 


2. Mathematical Model 


The starting point for the present work is our re- 
cent model 4,5 that generalizes classical models 6,7 of a 
reactive liquid/gas interface by replacing the simple 
assumption of a fixed normal propagation speed with 
a reaction/pyrolysis rate that is a function of the lo- 
cal pressure and temperature. This introduces im- 
portant new sensitivity parameters that couple the 
local burning rate with the pressure and tempera- 
ture fields. Thus, it is assumed, as in the classical 
models, that there is no distributed reaction in ei- 
ther the liquid or gas phases, but that there exists 
either a pyrolysis reaction or an exothermic decom- 
position at the liquid/gas interface that depends on 
local conditions there. In its most general form, the 
model includes full heat and momentum transport, 
allowing for viscous effects in both the liquid and 
gas phases, as well as effects due to gravity and sur- 
face tension. For additional simplicity, however, it 
is assumed that within the liquid and gas phases 
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separately, the density, heat capacity, kinematic vis- 
cosity and thermal diffusivity are constants, with ap- 
propriate jumps in these quantities across the phase 
boundary. 

The nondimensional location of this interface 
as a function of space and time is denoted by x 3 = 
$ s (zi,x 2i £), where the adopted coordinate system 
is fixed with respect to the stationary liquid at x 3 = 
-oo (Figure 1). Then, in the moving coordinate 
system x = x u y = x 2 , z = x 3 - $,(xi,x 2 ,t), in 
terms of which the liquid/gas interface always lies 
at z = 0, the complete formulation of the problem is 
given as follows. Conservation of mass, energy and 
momentum within each phase imply 


V • v = 0 , 0, (1) 


9 © 

_ 9£, 

9© 

dt ' 

dt 

dz 

9v 


dv 

dt ' 

dt 



+ v • V9 = {*}v 2 ©, z< 0,(2) 
+ (v • V)v = (0, 0, — FY _1 ) 


-{p-.} Vp+ {£r'J VV ' *5°' 


( 3 ) 


where v, © and p denote velocity, temperature and 
pressure, respectively, Pr t and Pr g denote the liquid 
and gas-phase Prandtl numbers, p, X and c (used be- 
low) are the gas-to-liquid density, thermal diffusivity 
and heat-capacity ratios, and Fr is the Froude num- 
ber. Here, the nondimensional variables have been 
defined in terms of their dimensional counterparts 
(denoted by tildes) as 


v 

v = — 

U 


P = 


$ a = -4 




Xi 


PiU 2 ’ 

t - — 
~ A, 


© = 


T-r u 

f a -T u 


XiU 

Xi ~ A , 


( 4 ) 


where U is the reference propagation speed of the 
interface for the case of steady, planar deflagration. 
The nondimensional parameters, some of which first 
appear below in the conditions at the gas/liquid in- 
terface (at z = 0), have been defined as 


liquid and gas-phase kinematic viscosities. We note 
that the Froude number Fr represents the ratio of 
inertial to buoyant (gravitational) forces, and, for fu- 
ture reference, that pXPr g = pPr t , where p = p g jpi 
is the gas-to-liquid viscosity ratio. In addition, we 
introduce the nondimensional mass burning rate 

^(piz=o + ,©u=o) = , (6) 

which is assumed to depend on the local pressure and 
temperature fields in the vicinity of the liquid/gas 
interface. Clearly, A = 1 for the case of steady, pla- 
nar burning, but perturbations in pressure and/or 
temperature result in corresponding perturbations 
in the local mass burning rate. 

Equations (1) - (3) are subject to the boundary 
conditions 

v = 0 , Q = 0 at z = -oo , 

© = 1 at z = +oo , ©| z=0 - = ©| z=0 + ^ ‘ ' 

and appropriate jump and continuity conditions at 
the liquid/gas interface. The latter consist of con- 
tinuity of the transverse velocity components (no- 
slip), 

n s x v_ = n s x v + , (8) 

where v± ~ v| z _ 0 ±, conservation of (normal) mass 
flux, 

n„-(v_ - pv + ) = (1 - p)S(&,)— - , (9) 

the mass burning rate (pyrolysis) law, 

9$. 

n s -v_ - S(9.)— = A(p\ z=0+ , 0| z=o ) , (10) 

conservation of flux of the normal and transverse 
components of momentum, 

p \ z = o ~ ~p|z=o+ = n s -(v_ - pv + )S($ 5 ) — -£ 

at 

+ [P v +(n s • v+) - v_(n s • v_) 

-pXPr g e+ • n a + Prie_ • n„] 


■b 

1) 

5 '| 5 " 

|*Eh 

II 

b 

eo Li* 

I 1C?) 

II 

II 

c-4. 

Cl 

II 

< 

Pr '=l' 

Pr - ^ 

r ’ 9 r y 

A 9 


- ( 5 ) 
where T u is the unburned (liquid) temperature at 

z = — oo and T a is the adiabatic burned (gas) tem- 
perature at z = -boo, 7 is the coefficient of surface 
tension for the liquid surface, and vi and v 9 are the 



I r, , /^M 2 ! _ 9% \ 

dy 2 [ \ 8x ) \ dx dy dxdy J 1 

r 

n s x [pv + (n. • v + ) - v_ (n s • v_ ) 

+ (v_ -pv + )5($,)^i] = (12) 

n s x (pAPr p e_|_ • n s - Pr z e_ • n s ) , 
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and conservation of heat flux 
n s • (cpAV0| 2=o + - V0| 2=o -) = 

n s -[(cpv+ - v_) 0 | z=o + c(a u pv + - v_)] 

+ [(1 - cp)©U=o + c(l - , 

(13) 

where c = c/(l - cr„), e is the rate-of-strain tensor 
(e± =e| 2==0 ±), 7 is the surface-tension coefficient, 
a u is the unburned-to-burned temperature ratio, and 
i?(^ 5 ) and the unit normal n s are defined as 

S(<f>,) = [1 + {d$ 3 /dx) 2 + {dsjdy ) 2 ] " 1/2 , 
n s = (~d*,/dx, -9$ s /9y, 1)S(*.) . 


respectively. Finally, we note the influence of viscos- 
ity on the jump in pressure across the liquid/gas in- 
terface through the appearance of the Prandtl num- 
bers Pri and Pr g in Eq. (11). We also remark that 
the momentum equation (3) suggests, as noted by a 
reviewer, a singularity in the pressure distribution at 
the interface itself due to the discontinuity in veloc- 
ity there. However, this singularity does not affect 
the jump in p across z = 0 that follows from conser- 
vation of the normal component of momentum flux, 
Eq. (11), which is all that is required to close the 
present model. 

A nontrivial basic solution to the above prob- 
lem, corresponding to the special case of a steady, 
planar deflagration, is given by 


Here, the factor multiplying 7 in Eq. ( 11 ) is the cur- 
vature —V • n of the liquid/gas interface in the mov- 
ing coordinate system, and the corresponding ex- 
pressions for the gradient operator V and the Lapla- 
cian V 2 in this system are given by 


/ d d$ a d d d$ s d_ d \ 

\ dx dx dz 1 dy dy dz' dz) 


(15) 


and 


V 2 


d 2 d 2 

dx 2 + dy 2 


1 + 


f d$g 

V dx 


) s +(t) 


d 2 

dz 2 


gg, d 2 99, d 2 

dx dxdz dy dydz 

_ d 2 $ 3 \ d 

\ dx 2 dy 2 / dz 

(16) 

However, the vector v still denotes the velocity with 
respect to the (x\ 1 X 2 ,xs) coordinate system. We re- 
mark that the mass burning rate A may be typically 
decomposed as 


^(p|z-o+»©U-o) = 


^Cp|z=o+,Q|z=o)exp 


' N{ 1-0(61^0-1) 1 
cr u + (1 - cr lt ) 0 | 2= o 


U?) 

where i is a rate coefficient and N = E/R°T a is 
the nondimensional activation energy (E), However, 
this more explicit representation will not be needed 
in the stability analysis that follows. Instead, re- 
sults may be expressed in terms of the pressure and 
temperature sensitivities, defined as 


A P = dA/dp\ e=l p=0 , 

A© = 9A/3©| e _ l p _o = Af(l — cr u ) + A© , 


= 0 °(z) = | 1 ] 

v°-(0 1 0,w°) 1 w° = | 



— Ft 1 z + p 1 
—pFr~ 1 z, 


, z < 0 
z > 0 , 

0 , z < 0 

p~ l — 1 , z > 0 , 

— 1 , z < 0 
z > 0 . 


(19) 


The linear stability analysis of this solution now pro- 
ceeds in a standard fashion. However, owing to the 
significant number of parameters, a complete analy- 
sis of the resulting dispersion relation is quite com- 
plex. Realistic limits that may be exploited to facil- 
itate a perturbation analysis of the dispersion rela- 
tion include p < 1 , \i 1 , and in the microgravity 
regime, Fr~ l 1 . In contrast, the earlier classi- 
cal studies considered special limiting cases and/or 
assumptions. Thus, in the study due to Landau , 6 
viscosity was neglected and the effects of gravity (as- 
sumed to act normal to the undisturbed planar in- 
terface in the direction of the unburned liquid) and 
surface tension were shown to be stabilizing, leading 
to a criterion for the absolute stability for steady, 
planar deflagration of the form (in our nondimen- 
sional notation) 47 Fr -1 p 1 /(l — p) > 1. In the study 
due to Levich , 7 surface tension was neglected, but 
the effects due to the viscosity of the liquid were 
included, leading to the absolute stability criterion 
Fr' x Pri{Zp)^ 2 > 1. Thus, these two studies, both 
of which assumed a constant normal burning rate 
(A = 1), demonstrated that sufficiently large values 
of either viscosity or surface tension, when coupled 
with the effects due to gravity, may render steady, 
planar deflagration stable to hydrodynamic distur- 
bances. In the present work, these results will be 
synthesized and extended to the more realistic case 
of a nonconstant burning rate (Le., A p ^ 0) in both 
normal and reduced gravity regimes. 
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3. The Hydrodynamic Linear Stability Problem 


With respect to the basic solution (19), the var- 
ious perturbation quantities <j> 3 (x,y,t), n(x,y,z,t), 
C(x,y,z,t) and 6(x,y,z,t) are defined as 


$ *=$2(0 + <ftn v = v°(z) + u , 

p=p°(z) + C. 0 = e°(z) + 9 + 4> s ^-. 


(20) 


Substituting Eqs. (20) into the nonlinear model (1) 
- (16) and linearizing about the basic solution (19), 
we obtain a problem for the perturbation quantities 
as 


du i du2 du3 

dx dy dz 


z ^0, 


( 21 ) 


r 1 1 du du 
t pJdt dz 



dx’dy 



d<Ps dC\ 
dy ’ dz) 


( Pri wd 2 u <9 2 u 9 2 u\ 

+ \ p\Pr g J + + lb 2 ) ’ 

\p)dt dz~\ o J 

+ flW*l M M) 

\pXi\dx 2 dy 2 dz 2 )’ > 


z > 0, 
(22) 

(23) 

0, 


u = 0 , ^ = 0 at 2 = — oo , 

6 = 0 at z = -f oo , 8\ z =o+ ~ = 4>s , 


p\Pr 9 ( 
- Pri 


' du 2 

, & z 
/ du 2 
\!h 


z=o+ d V 
du 3 

z=o- dy 

d<t> 


z= 0+ 


z=0- 


( 31 ) 


f + “ 2 U + - 2 U-. 


cpX 


86 

dz 


r=0+ 


86 

dz 


” ^[2=0+ + 8\z=§~ ~ 

z= 0- 


cpu 3 \ z=0+ -(l+c)u 3 | z=0 . +[l+c(l-p)]^, 

(32) 

where c = c(l - cr u ) \ 

In the present work, we shall focus, using our ex- 
tended model described above, primarily on hydro- 
dynamic (Landau) instability. Thus, in the linear 
stability analysis, we retain only the pressure sen- 
sitivity Ap in the perturbation of the pyrolysis law 
(10), neglecting the temperature sensitivity Aq. The 
latter assumption thus filters out reactive/diffusive 
instabilities associated with the thermal coupling of 
the temperature field, 4 ’ 5 but greatly facilitates the 
analysis of instability due to hydrodynamic effects 
alone. We note that the mass burning rate of many 
propellants has been shown empirically to correlate 
well with pressure. 

Nontrivial harmonic solutions for u and (, 
proportional to e tUii+ikl x + ik *y } that satisfy Eqs. (21) 
- (22) and the boundary/boundedness conditions at 
z = ±oo are given by 


«i| z=0 - -^i| z=0+ =(p" 1 -l)^, (25) 

^U_-u 2 | z=0+ =(p- 1 -l)^ > (26) 

tt3U 0 - -ptx 3 | z=0+ =(1 (27) 

i dd> s . , 

U3 lz=o- gf ~ ^pClz=o+ + j4e0| z=o + , (28) 


C| z =o- - Clz=o+ = (2 - p)u 3 \ z=Q+ - U3| z=0 






z=0 


z=0+ 


d<p 3 


(•j ^(^ 2< ^ S , d 2( t > s\ 

plPr,( a i i| + £i ) 

\ lz=0+ z=0+/ 

-Pn( a 4 i| + 


V dz lz= 0 - 


z=0- 


d<t>. 


(p_1 _1) '^ + Ul lz=o + ~ U A 


z— 0“ 


(29) 


(30) 


= giut+ikiX+ikiy 

(33) 

i iut+ikix+ik 2 y / b\e kz — Ft ^ 

(34) 

\ b 2 e~ kz - pFr ~ l , 


u e iut+ikix+ik 2 y f b 2 e qz — iki(itu 4- k) l b\e kz 

\ b\e TZ — ik\(iujp — k)~ l b 2 e ~ kz , 

(35) 

u giut+ikix+ikzy f b$e qz — ik 2 (iu) 4 k)~ l bie kz 

\ bee rz — ik 2 (iup - k)~ l b 2 e ~ kz , 

(36) 

_ giut+ikix+ikiy ( bje qz — k(iu 4* k)~*b\e kz 

1 b$e rz -f k(iujp - k)'~ 1 b 2 e~ kz 1 

(37) 

for z > 0, where we have normalized the above so- 
lution by setting the coefficient of the harmonic de- 
pendence of 4> 3 to unity. Here, the signs of 1 fc x and 
k 2 may be either positive or negative and we have 
employed the definition k = + Ar|, and q and r 

are defined as 
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2 Pri q = l + y/l + 4Pr((zw + Pr, it 2 ) , (38) 


2 fiPri r = 1 - y/l + 4p,Pri(iup + pPri k 2 ) , (39) 

where we have used the fact, noted below Eq. (5), 
that pXPr g = jiPri. 

Substituting this solution into the interface con- 
ditions (25) - (31) and using Eq. (21) for z > 0 yields 
nine conditions for the eight coefficients b\ - b% and 
the complex frequency (dispersion relation) iu(k). 
In particular, these conditions are given by 


ikibs 4- ik 2 bs 4- qb 7 = 0 , (40) 


ifci& 4 4- ik 2 b 6 + rb s = 0 , (41) 


63 - 

iki 

- . . , 0 1 - 

- 64 4- 

ik\ 

62 = 

(l-i 

[)iki , 

(42) 

xu 4* k 

itup — k 




b 5 ~ 

ik 2 
, , 0\ 

- 4- 

ik 2 

b 2 — 

(1-: 

[)ik 2 , 

(43) 


xu 4- k 

iup — k 




b 7 - 

k h 

- . , , - 

-pb 8 

pk 

— b 2 z 

= (i- 

p)iu, 

(44) 


VjJ 4- k 

iup — 
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b 7 - . ^ -bi - A v b 2 = iu - pFr l A p , (45) 

to; 4- k 


[ 1 + s4e<**<-» 
-[ 


iu) + k 

1 + 


6 i 


-(2 kpPri +2 — p) 


b 2 


iuip — A: 

+ (1 - 2Pn g)fe 7 - (2 - p- 2 pPn r)b s = 

(1 - p)(Fr ~ 1 - iu>) + 7 fc 2 , 


(46) 


(/iPri r - 1)64 + (2 kpPri + 1) . — , b 2 

XU p K 

+ ikipPri b a + (1 - Pr t q)b 3 - ikyPri b 7 (47) 
+ (2Pr t k - 1) - ~-r h = (1 - , 

{pPri r - 1 )& 6 + (2 kpPri 4- 1 

ILjp — K 

4- ik 2 pPri bs + (1 - Pri q)b 5 - ik 2 Pri b 7 (48) 

+ (2 Pr, k - l)-r— ^-7-61 = (- - l)ik 2 . 

VjJ + k \p J 


While the above problem is linear in the coefficients 
61 - 6g, the relationship for iu is highly nonlinear. 
Accordingly, we ultimately will exploit the fact that 
the gas-to-liquid density and viscosity ratios p and p 
are typically small, as is Fr*" 1 in the case of reduced 
gravity, and seek an asymptotic representation of the 
dispersion relation in this parameter regime. As we 
shall see, this distinguished limit simplifies the prob- 
lem by inducing three distinct wavenumber regimes 
that can be analyzed individually. Consequently, it 
is possible to develop separate perturbation expan- 
sions for the neutral stability boundaries in each of 


these wavenumber regimes, which can then be com- 
bined using asymptotic matching principles to pro- 
duce a uniformly valid composite expression for each 
boundary. 



Hydrodynamic neutral stability boundaries in the 
limit of zero viscosity. 


4. Hydrodynamic Instability: Zero-Viscosity Limit 


Although both liquid and gas-phase viscous ef- 
fects will be shown to be comparable in general to 
those due to surface tension, the inviscid limit pro- 
vides a particularly tractable limiting case and moti- 
vates the scalings to be used in the asymptotic anal- 
ysis of the fully viscous problem. Thus, in the limit 
of zero viscosities (Pri = Pr g = 0), our extended 
model differs from the classical one due to Landau 6 
only in the local pressure sensitivity of the normal 
burning rate. In that limit, the dispersion relation 
with respect to infinitesimal hydrodynamic distur- 
bances proportional to e lwt±lk x > where k and x are 
the transverse wavenumber and coordinate vectors, 
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respectively, is determined as 5 

A p p 2 (l - P)(iu) 3 ~ P 2 k [1 + p + A p ( 1 - p)} (iu) 2 
+ P k {A v [p(l + p)Fr~ l + pk 2 7 - (1 - p)k ] 

-2 pk} (iu>) 

+ A p fc 2 [(l - p)(2 - p)k + p 2 { 3 - p)Fr~ 1 + p 2 k 2 1 \ 

+ pk 2 [(1 - p)k - p{ 1 - p)Fr~ l - pk 2 7] = 0 . 

( 49 ) 

The corresponding neutral stability boundaries are 
easily determined explicitly and are exhibited in Fig- 
ure 2. Steady, planar burning is always unstable for 
positive values of A p , but in the region A p < 0, there 
exist both cellular (u; = 0) and pulsating (a? ^ 0) 
stability boundaries A p (k\ p, 7, Fr” 1 ) given by 

a = p M 1 ~ ^ Fr ~ 1 + Z t 1 ~ < n 

p p 2 (3-p)Fr~ 1 + p 2 ^k 2 + (1 -p)(2 -p)k ~ 
oj = 0 , 

(50) 

and 


respectively, where k = |k|. Steady, planar com- 
bustion is therefore stable in the region A p < 0 
that lies between these two curves, where the pul- 
sating stability boundary is a straight line in the 
( A p , k) plane and the cellular stability boundary is 
a curve which lies above the former in the region 
A p > — p/(2 — p). The shape of the latter bound- 
ary depends on whether or not the parameters Fr” 1 
and/or 7 are zero. In the limit that the product 
7 Fr” 1 approaches the value (1 — p)/4p 2 from be- 
low, the cellular stability boundary recedes from the 
region A p < 0. For 7 Fr” 1 > (1 — p)/4p 2 , the sta- 
ble region is the strip -p/( 1 — p) < A p < 0. Thus, 
when A p = 0, the classical Landau result for cellu- 
lar instability is recovered. However, even a small 
positive value of A p renders steady, planar burning 
intrinsically unstable for all disturbance wavenum- 
bers, regardless of the stabilizing effects of gravity 
and surface tension, since there always exists a pos- 
itive real root of Eq. (49) for A v > 0 (for small p, 
this root is given by iu ~ k/A p ). This result may 
be anticipated from quasi-steady physical consider- 
ations. That is, a burning velocity that increases 
with increasing pressure is a hydrodynamically un- 
stable situation, since an increase in the burning ve- 
locity results in an increase in the pressure jump 
across the liquid/gas interface, and vice-versa. How- 
ever, a sufficiently large negative value of A v re- 
sults in a pulsating hydrodynamic instability, the 


existence of which was a new prediction for liquid- 
propellant combustion. Zero and negative values of 
A p over certain pressure ranges are characteristic of 
the so-called “plateau” and “mesa” types of solid 
propellants, 10 as well as for the HAN-based liquid 
propellants mentioned above. 3 

Of particular interest in the present work is 
the cellular hydrodynamic, or Landau, instability of 
liquid-propellant combustion in the limit of small 
gravitational effects (in microgravity, for example). 
In this limit, the shape of the upper hydrodynamic 
stability boundary in Figure 2, corresponding to the 
classical type of instability that was first described 
by Landau, 6 clearly approximates the Fr” 1 = 0 
curve except for small wavenumbers, where, unless 
the inverse FYoude number is identically zero, the 
neutral stability boundary must turn and intersect 
the horizontal axis. Consequently, the neutral sta- 
bility boundary has a minimum for some small value 
of the transverse wavenumber k of the disturbance, 
implying loss of stability of the basic solution to long 
wavelength perturbations as the pressure sensitivity 
A p defined above decreases in magnitude. This, in 
turn, suggests a small wavenumber nonlinear stabil- 
ity analysis in the unstable regime, which generally 
leads to simplified nonlinear evolution equations of 
the Kuramoto-Sivashinsky type for the finite ampli- 
tude perturbations. 11,12 

To establish the nature of hydrodynamic insta- 
bility in the microgravity regime in a formal sense, as 
well as to introduce a formalism whereby the anal- 
ysis of the fully viscous case may be simplified, we 
may realistically consider the parameter regime p 
1, Fr -1 1, with Fr” 1 ~ p. For example, typical 
values are p ~ 10” 3 — 10” 4 , liquid thermal diffusiv- 
ity A 1 -0.1 m 2 /sec, and the steady, planar burning 
rate U — 1 — 10 cm/sec depending on pressure. 3 
Hence, from the definition Fr” 1 = gXi/U 3 y we con- 
clude that Fr” 1 ~ p implies that the dimensional 
gravitational acceleration g < 10” 5 m/sec 2 , which, 
roughly speaking, marks the onset of the micrograv- 
ity regime. Thus, introducing the bookkeeping pa- 
rameter e « 1, we define scaled parameters p*, p* 
and A* according to 

P = P*t, A p = A ;e , Fr~' = { 9 g . € , (52) 

with 7 remaining an 0(1) parameter. Here, the 
lower scaling on Fr” 1 corresponds to the reduced 
gravity limit, whereas the upper definition indicates 
the normal gravity case. In both regimes, it is read- 
ily seen from Eq. (50) that there are three distinct 
wavenumber scales that must be considered in order 


k 2 k 

7 4 - - , 

1 -p pj 

( 51 ) 


A p = , 

1-P 


= ‘[tt 


-Fr' 1 + 


1 


to properly account for the relative magnitude of the 
various terms in the numerator and denominator in 
the expression for A p . In particular, in addition to 
the 0(1), or outer, scale k, there is an inner scale fc* 
and a far outer scale k / that are respectively defined 
as 


/ k/e, Fr~ l -0(1) 
{ k/e 2 , Fr~ l ~0{e) 


k f = ke. (53) 


In the thin inner and thick far outer regions, we thus 
readily obtain 


A 


A P 




/ P*{P*9 ~ ki)/2ki 


(54) 


and 

a;~a;v ] ( 55 ) 

respectively. Each of these two expansions may be 
matched to the 0(1) outer expansion, which is given 

by 



and thus a uniformly valid composite expansion, de- 
noted by j4*^(/c), may be constructed as 


a; ( c) ~ 

a;" 


+ A* (o) + Al u > - lim Al w - lim A 


«-(/) 


(0 


*■(/) 


ki- 


kf—+ 0 




1 , 1 ,2 , 

~2 P + 2 ep lk + 


I ep^g/M 

l «V V/2* 


(57a) 

where the definitions of ki and k j have been used to 
express the final result in terms of k (Figure 3). In 
terms of the original unsealed parameters, Eq. (57a) 
beomes 


A p 




P p^yk p^Fr 1 

2 2 2k 


(57 b) 


which eliminates the small scaling parameter e from 
the composite expansion. Noting the simplicity of 
the asymptotic result (57) relative to the exact ex- 
pression given in Eq. (50), it is easily seen that the 
hydrodynamic stability boundary in the regime con- 
sidered here lies in the region A* < 0, intersecting 
the A* = 0 axis at 


and at 


k = ho — 


~ ~ — » 1 
P* 

(58) 

^ 2 «i, 
P 9 t 

(59) 


with a single local minimum at 


corresponding to the critical value 




ep* 2 \/57 

ep' 2 y/eg*j 


(61a) 


Thus, in terms of the unsealed parameters, instabil- 
ity first occurs for disturbances at the critical value 

A p = A p + p 2 VlFr - 1 . (616) 


The essential difference, as illustrated by Figure 
3, between the normal and reduced gravity limits in 



Asymptotic representation of the cellular hydro - 
dynamic neutral stability boundary in the limit of 
zero viscosity. The upper (lower) solid curves cor- 
respond to the two cases described by Eqs. (57) for 
normal and reduced- gravity, respectively (curves are 
drawn for the case e = 0.04, p* = 1.0, g = 2.5, 
9 ' = 10 ). 
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the realistic parameter regime considered here is that 
in the latter instance, gravity is only capable of sta- 
bilizing disturbances whose wavenumbers are very 
small, 0(e 2 ), whereas in the former case, gravity 
is sufficiently strong to stabilize disturbances whose 
wavenumbers are 0(e). As a consequence, hydro- 
dynamic instability becomes a long-wave instability 
phenomenon in the reduced gravity regime consid- 
ered here, since the most unstable wavenumbers are 
0(y/c) ) rather than 0(1), in that case. We note that 
the role of surface tension, on the other hand, serves 
to stabilize only short-wave disturbances, which are, 
to this order of approximation, unaffected by the 
magnitude of the gravitational acceleration. 


5. Hydrodynamic Instability: Fully Viscous Case 


Guided by these results for the inviscid case, the 
linear stability analysis may be extended to include 
the effects of viscosity as follows. Retaining the scal- 
ings (52), we have noted that p\Pr g — pPru where 
p = ii g / pi is the gas-to-liquid viscosity ratio. Thus, 
it is reasonable to treat Pri = P as an 0(1) param- 
eter, and to consider the limit 

p = p* e <C 1 . (62) 


Introducing the scalings (52) and (62) directly into 
Eqs. (40) - (48), approximate solutions for the coef- 
ficients may be sought in the form of appropriate ex- 
pansions in powers of e. As suggested by the inviscid 
analysis, it is necessary to construct separate expan- 
sions for the inner, outer and far outer wavenumber 
regimes defined above. Since we are interested in ob- 
taining the hydrodynamic cellular stability bound- 
ary, we may set iu equal to zero directly and seek 
solutions for b\ — bs and A*, where the latter ex- 
pressed as a function of k will describe the stabil- 
ity boundary. Forms for the various expansions are 
suggested in part by the corresponding solutions for 
these coefficients in the inviscid limit (for which ex- 
act solutions are analytically tractable 5 ). 

In the 0(1) wavenumber regime, the expansions 
for q and r are, from Eqs. (38) and (39), 


r ~ r, c + • • • 
q ~ qo H i 


n = n'Pk 2 , 

l + vT+Ip^F 


qo = 


2P 


(63) 

(64) 


and appropriate expansions for the coefficients 6, — 
i> 8 and A v are given by 


bi ~ V 1 + b| 0) + b^e + ■■■ , i = 2, 8(65) 
6i~&l° ) +6{ l) e + ■ , i = 1,3, 4, 5, 6, 7,(66) 

A p = A; (o) £ ~ e(A' 0 {o) + A; (o) £ + •••)• (67) 


Substituting these expansions and the scalings (52) 
and (62) into Eqs. (40) - (48), and collecting coef- 
ficients of like powers of e, we obtain a sequence of 


algebraic problems for the coefficients in Eqs. (65) 
- (67). In particular, from either of Eqs. (42) and 
(43), we obtain 

<4' X) = -fc/p* , (68) 

and the leading-order versions of Eqs. (44) - (46) 
are given by 

4 o) - 4 0) - P *4 _1) + p*4 _1) = o , (69) 

6 W_ 6 i° ) _ j 4 ,(o) 6 (- 1) =0 , (70) 

&2 _1) - 2f4 -1) = 0 . (71) 

Subtracting Eqs. (69) and (70), and using Eqs. (68) 
and (71) for 63^ and 6g -1 \ we obtain the result 

A'o {o) = ~\p* • ( 72 ) 


Equations (67) and (72) are equivalent to the invis- 
cid result (56) in the outer wavenumber regime, and 
thus, to the leading order of approximation, viscos- 
ity does not affect the neutral stability boundary for 
0(1) disturbance wavenumbers. 

A corresponding result is obtained in the inner 
wavenumber regime defined by the first of Eqs. (53). 
In particular, the expansions for q and r are now 
given by 


r ~ p 


'Pk 2 | 


€* + ■■■ 
£ 5 + • • • 




€ 2 + ' 
e 4 + ' 


(73) 

and the expansions for the coefficients 6, - bs and 
A p may be sought in the form 


bi 

bi 

A p 



b W + b? ) e+... 
6i 1) e + fcJV + ... 1 

i = 2 , 8 , 

6 u>’ e + 6 (V + ..- 

, i = 1, 3, 

6< 15 £ 2 + 6- 3) £ 3 + ••• 

1 * *1 

;W£~e(A5«+Ai ( ’ 



(74) 


4, 5, 6/75) 
(76) 


Substituting these expansions and the scalings (52) 
and (62) into Eqs. (40) - (48) as before, we again 
obtained a much-simplified set of algebraic equations 
at each order. In the reduced-gravity regime, for 
example, we obtain from Eq. (40) and either of Eqs. 
(42) or (43) that 


6 


( 2 ) 

7 


= 0 , 



(77) 
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while Eqs. (44) - (46) give 


b {2) 

°7 


b[ 2) - P*bi 1] + p*bi 1] = 0, 

(78) 

-b^-Afh^=0, 

(79) 

6^ 1) -26< 1) = g\ 

(80) 


h&) 


An identical set is obtained in the normal-gravity 
regime in terms of the leading-order coefficients in 
the expansions (74) and (75) for that case, with the 
exception that g* is replaced by its unsealed coun- 
terpart g in Eq. (80). Subtracting Eqs. (78) and 
(79) and using the results (77) and (80), we thus 
obtain 

\P m (p*9-ki)/2ki , gl x 

^ \ p*{p*g* — ki)/2ki ’ W 

which is again equivalent to the inviscid result (54). 
Thus, to leading order, neither the inner nor the 
outer wavenumber regimes are influenced by viscous 
effects. Indeed, to a first approximation, viscosity is 
only significant for large- wavenumber disturbances, 
as we shall now demonstrate. 

In the far outer wavenumber regime defined by 
the second of Eqs. (53), the expansions for q and r 
are given by 

1- (1 +4/i* 2 P 2 fc2) 1 /2 

r ~ r <- 1 > € r <-«“ 2 JFP ’ 

(82) 

<7(-i) = kf , (83) 

and appropriate expansions for the coefficients &i - 
b$ and A p are given by 

fc<0) 


i)C 1 + 


6 j ~6i“ I) e“ l + 6} 0, + ---, < = 1,3,5, 7, (84) 

bt ~ b[~ 2) e- 2 + fej" V 1 +■■■, i = 2, 4, 6, 8, (85) 
A p = A; (f) e ~ e(A* 0 U) + A\ U) t + •••)• (86) 

Substituting these latest expansions into Eqs. (40) 
- (48), we obtain the leading-order system 


ik lf b[ 


+ ik 2 j 65 "I" Q( — i) bj 0, 


ikif b[~ 2) + ik 2 f b ( ~ 2) +r/_,\ b[~ 2) = < 


-b{~ 2) - 


2f°6 + r (-l) u & 

ik i f . (_^\ ik 

I 


£!//>(— 2 ) 
kf 2 


1/ 


P* 


-b { ~ 2) - 

°6 


ik2f u (-2) _ ^2/ 

Do 


7 u 2 “ * » 

k f p 

°8 


b^-bi-V-p'b'fV+p'bi- 2 ^ 0 , 


(87) 

( 88 ) 

(89) 

(90) 

(91) 

(92) 


2 k f Pb[- 1] + (1 + 2k ft i*P)b { 2 - 2) - 2P< 7( _ 1) 6< l) 

-2(l-M*Pr ( _ 1) )6^ 2) =7fc? 

ik 


(93) 


(M*Pr ( _D - 1)64 2) - (2k }t x'P + 1 




+ ik lf p.*Pb , f 2) - Pq^bf 1 ’ + 2 ik lf Pb\ 




k(-D 




(M*Pr ( -i) - 1)^ _2) - (2fc/M*P + 1 )^^ 

Kf 


ik 


(94) 


Vt(- 2 ) 


+ ik 2f ^Pbi~ 2) - Pg ( _i)6r i; + 2ik 2 fPb\ 

(95) 

where, analogous to the scaling for k in the far outer 
regime, k\j = fcie and k 2 f = £2^ Equations (87) - 
(95) may be solved for the leading-order coefficients 
in Eqs. (84) - (86) such that, after some algebra 
(Appendix A), we obtain the result 


(-1) 




A*o U) = —r—\{kf + r { - 1)) 

r (-i) L 

(r ( _i) 4- fe/)(l + 2p*Pfc/) +k f 


P*7 r (-i) - 2 P*Pkj ~ 2(1 - £i*Pr-(_i))J ’ 

(96) 

which may be rewritten as 


a+U) 

^0 




4- 


2 p>*P [1 + fc/(p*7 + 2m*P + 2p*P)] 
4p*P(l + p*Pfc,) - [1 - P(fe/)] (p*7 + WP) ’ 

P(fc/) = l+4p* 2 P 2 ^J . 

(97) 

It is readily observed that both the liquid and the 
gas-phase viscosities (through the parameters P and 
p*P, respectively) enter into this expression, reflect- 
ing an equal influence of viscous and surface-tension 
effects on the neutral stability boundary in the large- 
wavenumber regime. The equal importance of gas- 
phase viscosity relative to that of the liquid phase 
stems from the fact that gas-phase disturbances are, 
according to Eqs. (84) and (85), larger in magni- 
tude than those in the liquid phase, such that a 
weak damping of a larger magnitude disturbance is 
of equal importance to an 0(1) damping of a smaller 
magnitude disturbance. In the limit P — + 0, the in- 
viscid expression (55) is recovered. It is easily shown 
that 

lim A*W — n * 

k™0 A ° - 2 P ’ 
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(98) 



so that the far outer solution can be matched to the 
outer solution (72). As kj increases, increases, 
intersecting the Aq ^ = 0 axis at the value kf = kf 
given by 


Figure 4 that increasing the values of any of the 
parameters P , p*P or 7 serves to shrink the size 
of the unstable domain through damping of short- 
wave perturbations. The non-negligible effects of 


p*7 i 1 p*7 + 2p*p) ’ ^ 

which agrees with the inviscid result for p* <£ 1 in 
the limit p*P — ► 0. As this boundary crosses the 
horizontal axis, it must eventually reach a maximum 
and decrease to A = 0, since lim^^eo = 0. 
However, as indicated earlier in the discussion of 
the inviscid case, there is another root iu) 0( £ - 2 ) 
that takes on positive values for A* > 0 (Appendix 
B). Consequently, only that portion of the hydrody- 
namic stability boundary that lies below the hori- 
zontal axis ( Aq ~ 0) is relevant to our discussion. 

A uniformly valid composite expansion span- 
ning all three wavenumber regimes described above 
may be constructed as in Eqs. (57), giving the result 


a * 


(c) 


— p* + 


p* 2 r e 9 

2k l«V 


2 p*j£P_ [1 4- My 7 + 2 n*P + 2£P) 1 
+ 4p*P(l + ekp'P) - (p* 7 + 2fi m P) [l - R(ek)] ’ 

(100a) 

or, reverting to the unsealed parameters A p , p , p 
and Fr~ l ) 


2 ppP [1 + k(py + 2 pP + 2 pP)} 

+ 4/zP(l + kpP) “ (py + 2/zP)[l - P(*)] ’ 

(1006) 

which again eliminates the small bookkeeping pa- 
rameter £ from the problem. Both the normal and 
reduced gravity boundaries A p ^ c \k) are graphically 
exhibited in Figure 4 for various zero and represen- 
tative nonzero values of /P, P, Fr~ x and 7. It is 
readily seen that the essential qualitative difference 
between the normal and reduced- gravity curves is 
the location of the critical wavenumber for instabil- 
ity. Specifically, the minimum in the neutral sta- 
bility boundaries occurs for 0(1) values of k under 
normal gravity, and at k ~ 0(e l ^ 2 ) in the reduced- 
gravity limit considered here, as in the. inviscid case 
described above. Indeed, it may be shown that Eqs. 
(100) collapse to Eqs. (57) in the limit of zero viscos- 
ity (P = 0), but we again note that viscous effects 
in both the liquid (P) and gas (/PP) are compara- 
ble to surface-tension effects (7) in damping large- 
wavenumber disturbances. In fact, it is clear from 



Asymptotic representation of the cellular hydro- 
dynamic neutral stability boundary for the viscous 
case . The upper and lower sets of curves corre- 
spond to the normal and reduced- gravity regimes , re- 
spectively } in the asymptotic limit considered in this 
work (curves drawn for the case e = .04, p* = 1.0, 
g = 6.0, g * = 2.0/ The solid curves correspond to 
the inviscid limit (P = 0) with nonzero surface ten- 
sion (y = 2.5/ The dash-dot curves correspond to 
nonzero surface tension (y = 2.5) and liquid viscos- 
ity (P = 1.0/ but zero gas-phase viscosity (p*P = 
0/ The dash-dot-dot curves differ from the dash-dot 
curves by the addition of gas-phase viscosity (p*P = 
1.0/ and are similar to the dash- dot- dot- dot curves , 
where the latter correspond to larger viscosities (P — 
p*P = 2.0/ The dash- dot- dot- dot- dot curves corre- 
spond to a viscous case (P = /PP = 1.0/ but with 
zero surface tension (7 = 0/ so that, from Eq. (99), 
the curves do not intercept the A* = 0 axis. 
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gas-phase viscosity represents an important correc- 
tion to Levich’s original treatment 7 in which these 
effects were simply assumed to be small. The result 
(100a, b) thus synthesizes and significantly extends 
the classical Landau/Levich results, 6,7 not only in al- 
lowing for a dynamic dependence of the burning rate 
on local conditions in the vicinity of the liquid/gas 
interface, but also in its formal treatment of those 
processes (surface tension, liquid and gas- phase vis- 
cosity) that affect damping of large- wavenumber dis- 
turbances. 

6. Asymptotic Analysis of the Pulsating 
Stability Boundary 


As indicated previously, the existence of a non- 
stationary pressure dependence on the burning rate 
(i.e., A p 0) leads to the prediction of a pulsat- 
ing hydrodynamic stability boundary that is absent 
when such a pressure coupling is neglected, as in 
the original Landau/Levich theories. In the invis- 
cid case, this boundary (39) is a straight line that 
lies below the cellular boundary discussed above, but 
this is modified under the influence of viscosity, as 
we shall demonstrate. 

For the scalings (52) and (62) adopted in the 
preceding sections [in particular, for P ~ 0(1), p ~ 
0(e)], it turns out that, unlike the cellular stability 
boundary where viscous effects only have a leading- 
order effect in the far outer wavenumber regime, the 
effects of viscosity have a leading-order effect on the 
pulsating boundary for 0(1) wavenumbers as well. 
Thus, in the outer wavenumber region, we seek a 
solution for the dispersion relation in the form 

iuj ~ e (iu>o “1“ *F -h • * * ) , (101) 

where the leading-order term is suggested by the ex- 
plicit results for the inviscid case, 4,5 and the expan- 
sion in powers of e 1//4 is suggested by the expansions 
for r and q, which, from Eqs. (38) and (39), have 
the form 

r ~ r (1/ 2 )€ 1/2 + r {3/4) e 3/4 + r,e + • • • , 

f(i/ 2 ) = ~iw 0 p* , r (3 / 4 ) = -iu>ip’ , (102) 

T X = -ioJ 2 p* - (p’Pk) 2 , 


q ~ «?(_ i/4)C x ^ 4 + <7oe° H , 

r — 1 4 - iu>i / \/iujoJP 

Q(- 1 / 4 ) = V luJ o /P , Qo = 


(103) 


Corresponding expansions for the coefficients bi in 
Eqs. (40) - (48) are determined as 


6t=6 (-V 1 +6<- 3/4) e- 3 / 4 +6 l ( - 1/2) e - 1 / 2 
i = 1, 2, 8, 

6 i = 6 |- 1 / V 1 / 2 + 6 |- 1 / 4 ) e - 1 / 4 + --- 1 
i = 3, 4, 5, 6 , 


+ ’ ’ • I 

(104) 

(105) 


^ = 6H /4) e' 1/4 + 6fV + --- , i = 7, (106) 

where the leading terms in the expansions for &i, 
6 2l 64, be and b$ are consistent with the inviscid 
results 5 and the remaining coefficients appear only 
for nonzero values of P and are conservatively postu- 
lated to have the indicated expansions. Substituting 
these expansions into Eqs. (40) - (48) and equating 
coefficients of like powers of e, we obtain the leading- 
order equations/results 


ife 1 4- 1/2) +i* a 6i -1/2) +?(-i/4)4" 1/4) =°> ( 107 ) 



(108) 

(109) 

( 110 ) 

( 111 ) 

( 112 ) 


where Eq. (Ill) was obtained from the leading-order 
difference of Eqs. (44) and (45), and the remainder 
of the leading-order versions of Eqs. (40) - (48) give 
redundant results. Substituting Eqs. (109) - (111) 
into Eq. (112), we obtain 



(113) 


and thus (iu o) 2 < 0 for A* < ~p*/ 2, which es- 
sentially recovers the leading-order cellular stability 
boundary (56) for 0(1) wavenumbers, but gives no 
information on the pulsating boundary since iuo Is 
purely imaginary for Ap < —p*/ 2. Hence, stability 
in the latter region is determined by higher-order 
coefficients in the expansion (101) for iu>. 

Continuing with the analysis of the expanded 
forms of Eqs. (40) - (48), we obtain the second- 
order equations/results 

ik\ +<?(_!/ 4 ) b^ +<7o 4' 1/4) “ 0 , 

(114) 
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iki 64 1/4) + ik 2 b { 6 1/4) + rf i /2 ) b\ 


(- 3 / 4 ) 

’8 


+ r ( 3 / 4 ) &8 ^ = 0 ) 


6<- 3/4) = 0, 


kb[ 1/2) -g ( _i/4)4 V4) = (l - 


6 (-V4)_ T fc_ 6 (-3/4) =ict;i 
1 lUJo 

bt 3/4) = 0 , 

b[- 3/4) = 2b ( 8 - 3/4 >=0, 


6<- 1 / 2 )=6^- 1/2) = 


= 0, 


(115) 

(116) 

(117) 

(118) 

(119) 

(120) 
(121) 


where Eq. (117) was obtained from the sum of Eq. 
(42) multiplied by ik\ and Eq. (43) multiplied by z&2 
and the use of Eq. (107), Eq. (119) was obtained 
from the difference of Eqs. (44) and (45), and Eqs. 
(121) follow from Eqs. (47) and (48) in conjunction 
with Eq. (116). From these results and Eq. (107), 
we thus conclude that 

4" 1/4) = iuJ i = 0 , 6^" 1/2) = iu > 0 (l - , 

(122) 

where the fact that iuj\ = 0 implies the need to 
calculate UJ 2 to determine stability in the region 
A* < —p*/2. Proceeding in this fashion, we ob- 
tain from the next-order versions of Eqs. (44), (46), 
the difference of Eqs. (44) and (45), and the sum of 
Eq. (47) multiplied by ik\ and Eq. (48) multiplied 
by Z&25 the relations 


6< 0) - ^6) 
iu 0 


k fc(~l/ 2 ) , iuJ 2 + fc..(-l) 


+ 


(*w o ) 2 


kb i" 


p*b { s 4) + p*b ( 2 = iu > 2 , 


(123) 


k 
ted 0 


b[~ 1/2) + —{2Pk - 1)6) 


<-»+b<r 1/2) 


+ 


^ 0gl fe (-l)_ 2b (-l/2) = _ 


0 1 


(124) 


- p- b <f'/» + p- 6 <-™ + 


lUjQP ,(-l) 


+ - 4*62 1 ^ 2 ^ = —icoop* 


— iuop*b 8 ^ + kb 2 1/l2 ^ + 10/0/0*62 ^ 


+ iwob^ - (2P/c - 1) 


k 2 


tw 0 


b[~ 1] = 0, 


(125) 


(126) 


which, in combination with the expressions for b\ 
6<- 1) 


(-1) 


63 62 1//2 ^ and zcjq given above, constitute 


four equations for the four unknowns b[ 1/2 \ b^\ 
and iu 2- Solving these simultaneous equa- 
tions, we thus obtain 



iu 2 = A; 



(128) 


Equation (128) is the desired result, from which we 
conclude that iuj 2 > 0 for (A£///) 2 > 14-2P/c. Thus, 
in the region A* < 0, iu >2 vanishes on the boundary 


A* ~ —p*y/l + 2Pk , (129) 

which is a pulsating boundary (Fig. 5) since, from 
Eq. (113), iujQ is purely imaginary along this curve. 

Equation (129) is valid for 0(1) wavenumbers, 
but since it matches to the leading-order inviscid 
inner pulsating boundary A * = — p* as k — * 0, 
and becomes large in a negative sense as k becomes 
large, it is clear that Eq. (129) represents the pul- 
sating boundary for arbitrary wavenumbers. That 
is, for P ~ 0(1), the effects of (liquid) viscosity on 
the pulsating boundary are, to a first approxima- 
tion, absent for small wavenumbers, are first felt for 
0(1) wavenumber perturbations, and are sufficient 
to move this boundary to larger-magnitude values 
A* ~ 0(e^ 1 ^ 2 ) in the far outer wavenumber regime. 
In contrast, the cellular boundary (100) is unaffected 
for 0(1) and smaller wavenumbers, and is only modi- 
fied an 0(1) amount for 0(e _1 ) wavenumbers, where 
both gas and liquid viscosities are equally signifi- 
cant. Thus, the hydrodynamic pulsating boundary 
is more sensitive to liquid viscous effects than is the 
corresponding cellular stability boundary. In fact, 
even for smaller- magnitude viscosities such that P ~ 
0(e 1//2 (with p ~ O(e) as above), it is easily demon- 
strated (see Appendix C) that, in the A p regime con- 
sidered here, the pulsating boundary is eliminated in 
the far outer wavenumber region. However, for even 
smaller liquid viscosities such that P = Pc ~ O(e), 
it may be shown by a completely analogous calcu- 
lation (Appendix D) to the one given above that 
0(1) modifications to the pulsating boundary then 
occur in the far outer wavenumber regime accord- 

A 1 try 

ing to A* = -p*( 1 -{- 2 Pkf) , which, in terms of 

unsealed quantities, is the same as Eq. (129). 
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Figure 5 

Asymptotic representation of the pulsating hy- 
drodynamic neutral stability boundary for the viscous 
case (P > 0). The region between the pulsating and 
cellular boundaries (the latter are shown on an ex- 
panded scale in Fig. 4) is the stable region with re- 
spect to hydrodynamic instability. 

Other cellular and pulsating stability bound- 
aries are obtained for nonzero values of the temper- 
ature sensitivity parameter A©, and are thus of a 
reactive/diffusive nature since they arise from a cou- 
pling of the burning rate to the local temperature 
field. These have been analyzed in the realistic limit 
p 1 for the inviscid case, 5 and the generalization 
of these results to the fully viscous problem in both 
the normal and reduced-gravity regimes is currently 
under investigation. One important result obtained 
from the inviscid analysis is that a nonzero value of 
the thermal sensitivity A© appears to have little ef- 
fect on the hydrodynamic cellular stability boundary 
shown in Figures 2-4. This result suggests that the 
upper stability boundary in Figure 2, corresponding 


to the onset of steady cells on the propellant sur- 
face, is especially likely to be observable. Indeed, an 
analysis of nonlinear stability in the neighborhood 
of this boundary 8,9 not only confirms the existence 
of steady cellular structures above this boundary, 
but also demonstrates how the interaction of certain 
types of cellular modes can result in secondary and 
tertiary transitions to time-periodic motions 11,13,14 
that may also correspond to the sloshing type of be- 
havior observed in HAN/TEAN/water mixtures. 3 

7. Conclusion 

The present work has presented a formal asymp- 
totic analysis of hydrodynamic instability in liquid- 
propellant combustion for a surface model in which 
burning takes place at the liquid/gas interface. The 
model itself is based on a synthesized version of the 
classical models analyzed by Landau 6 and Levich, 7 
but generalized to allow a coupling of the mass burn- 
ing rate with the local pressure and/or temperature 
fields. 4,5 Focusing on the pressure-coupled version 
of the model, the realistic smallness of the gas-to- 
liquid density ratio was shown to be a convenient 
small parameter upon which to base an asymptotic 
treatment, resulting in three distinct wavenumber 
regimes such that different physical processes as- 
sume dominance in each. Both cellular and pul- 
sating hydrodynamic stability boundaries are pre- 
dicted, the former corresponding to Landau’s origi- 
nal notion of hydrodynamic instability, and the lat- 
ter representing a new prediction arising from the 
pressure dependence of the burning rate. For the 
cellular type of instability, it was shown that the 
gravitational acceleration (assumed to be normal to 
the undisturbed liquid/gas interface in the direc- 
tion of the liquid) is responsible for stabilizing long- 
wave disturbances, whereas surface tension and vis- 
cosity are effective in stabilizing short-wave pertur- 
bations. As a consequence, reduced gravity results 
in a shift in the minimum of the neutral stability 
boundary towards smaller wavenumbers, such that 
the onset of hydrodynamic instability, predicted to 
occur for sufficiently small negative values of the 
pressure- sensitivity coefficient A p , becomes a long- 
wave instability in that limit. An additional result is 
that gas-phase viscosity plays an equally significant 
role as does liquid viscosity in the large-wavenumber 
regime. This important effect, absent from previ- 
ous treatments, stems from the fact that gas-phase 
disturbances are larger in magnitude than those in 
the liquid phase. That is, although the gas-to-liquid 
viscosity ratio is small, a weak damping of a larger 
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magnitude disturbance is of equal importance to an 
0 ( 1 ) damping of a smaller magnitude disturbance. 
The inclusion of both viscous and surface-tension 
effects in a single analysis, which are of comparable 
importance for short-wave perturbations thus rep- 
resents an important synthesis of the classical Lan- 
dau/Levich theories of hydrodynamic instability. 

In the case of the pulsating hydrodynamic in- 
stability that arises from the pressure coupling that 
exists in the generalized model analyzed here, we 
find that neither gravity nor surface tension play a 
leading-order role, and viscous effects are the dom- 
inant stabilizing influence. Indeed, for 0(1) liquid 
Prandtl numbers, the stabilizing effects of (liquid) 
viscosity on pulsating instability are significant for 
disturbances whose wavenumbers are 0 ( 1 ) or larger. 
In contrast, it was shown that viscous effects are only 
significant for large- wavenumber disturbances in the 
case of cellular instability, where the influence of gas 
and liquid viscosity are comparable despite the small 
ratio of these two parameters. Although the onset 
of pulsating hydrodynamic instability is predicted 
to occur only for sufficiently negative values of the 
pressure-sensitivity coefficient A p) the persistence of 
the pulsating stability boundary (in the presence of 
viscous effects) for small wavenumbers suggests that 
it may be observable in those types of liquid propel- 
lants, such as those based on hydroxylammonium 
nitrate (HAN) and/or triethanolammonium nitrate 
(TEAN), that are characterized by negative pres- 
sure sensitivities over certain pressure ranges. How- 
ever, in attempting to interpret the sloshing behav- 
ior that has been observed during combustion of 
certain HAN/TEAN/water mixtures , 3 we note that 
nonsteady burning can arise via secondary and other 
higher-order bifurcations in the cellular instability 
region , 11,13,14 as well as from a primary crossing of 
the pulsating boundary. Hence, further measure- 
ments are generally needed to determine the precise 
origin of such behavior in any given experiment. 


Appendix A. Derivation of AJ^(A:/) 

The steps leading to Eq. (96) for the neutral 
stability boundary in the outer wavenumber regime 
are as follows. First, from Eqs. (89) and (90), we 
obtain expressions for 64 “ and b in terms of 
2 ^- Substituting these results into Eq. ( 88 ) then 


determines bg 2) in terms of 62 as well, so that 


&4~ 2) = 


n 

,(- 2 ) 

ikxjl- 

'2 

JL 

\ 

Kf 


,(- 2) 


*2 

kf 

k f 1 


r (- 1) ‘ 

2 

V 


(A.l) 


-b 


P * 


Using this last expression for bg ” 2 ' 1 in the equation 
obtained from subtracting Eq. (92) from Eq. (91), 
thus determines the neutral stability boundary in 


terms of b\ ^ as 


V r (-1) / r ( _ 1} 6< 2) 


(A.2) 


Using the last of Eqs. (A.l), we now solve Eq. (91) 
for by -1 ^ in terms of b[~ ^ and b 2 _2 \ 

4 - 1 ’ = A" (AL . + 1) - -3- , (A3) 

\ r (-l) / r (-l) 

and using this result along with the first two of Eqs. 
(A.l) in Eqs. (94) and (95) leads to the relations 


».C-a> 


L.(, + _t)_^.(A_ +2+ rh!>) 

V r (-i )J V r (-D k i )\ 

(-1) _ 


ikxf 3 +bl 

k } 

r (-l) 


k f p' k f \ 

L + b r (-D + ~ ’ 

r (-l) P \ r (-i )/ 

+ 4 ->- 

**(—&)■ 


(AA) 


r (-i) 


(A.5) 


Subtracting Eq. (A.5) from Eq. (A. 4), we obtain 




ih 


2/ 


b (_1) 
£3 

ikit 


which, when substituted into Eq. (87), gives 

4-" = ii-t-i = A,-. , 

7 ikif 3 ik 2 f 5 


(A. 6 ) 


(AS) 
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Finally, substituting Eq. (A. 3) and the last of Eqs. 

(A.l) into Eq. (93) and using the fact, from Eq. 

( — 2 ) 

(83), that q(-\) = A we completely determine b\ 
as 

(—2) _ P*7T~(-1) - 2p*Pfc; - 2(1 - M*^ >r (- 1) ) 

2 P * (»*(-i)+fc/)(l + 2p*Pfc/) + fc/ 

(>1.7) 

The remaining coefficients may now be solved for ex- 
plicitly. For example, b[~ 2 \ £>g and 63 are now 
determined according to Eqs. (A.l), while the use 
of Eq. (A. 7) for b and the first of Eqs. (A.6) for 
b?~ 1 ^ reduces Eqs. (A. 3) and (A. 4) to two simulta- 
neous equations for b[~ 1 '* and 63 l \ Furthermore, 
when Eq. (A. 7) is substituted into Eq.(A.2), we ob- 
tain the desired expression (96) for AjJ^. 


Appendix B. Intrinsic Instability for A* > 0 


A direct analysis of the dispersion relation (49) 
corresponding to the inviscid case indicates that, 
in addition to the roots that lead to the stability 
boundaries (50) and (51), there is an additional root 
that, for p ~ 0(e) < 1 and k ~ 0(1), is given by 
iu rsj e^fc/A*. 5 It may be shown that this result is 
valid in the inner and far outer wavenumber regimes 
as well, so that, at least in the inviscid case, the re- 
gion A* > 0 is intrinsically unstable. Based on the 
results in Section 5 for the viscous case, we may an- 
ticipate that the only possible modification to this 
result will occur in the far outer wavenumber regime 
k = fy/e, which we now consider. 

Thus, we seek a root of the dispersion rela- 
tion implied by the system (87) - (95) that has the 
asymptotic behavior in the far outer wavenumber 
regime given by 

iu ~ e 2 (ujq + *+■•••), (5.1) 

whereas, based on the inviscid results, the roots of 
the dispersion relation that leads to the cellular sta- 
bility boundary (97) in the far outer regime behave 
as 0(e -3 / 2 ). Thus, corresponding to Eq. (B.l), the 
expressions for r(_i) and Q(~i) in the expansions (82) 
- (83) for r and q are now given by 


whereas the appropriate expansions for the coeffi- 
cients bi are now written as 

fc i ~6|"V 2 + 6i‘V 1 + --- , t = 3,5,7, (5. 3) 

bi ~ 6|- 3) €" 3 + &! _2) e~ 2 + • • • , t = 1, 2, 4, 6, 8 . 

(5- 4 ) 

Hence, in place of Eqs. (87) - (95), the leading- 
order coefficients and the root iuj o are determined 
from Eqs. (40) - (48) by 


ikifb$ ^+^2/65 2) + <?(_!) &7 ^ — 0, {B-5) 

ik\f 64 + ik2f 4- ^(-1) ( B.6 ) 

— n (B.7) 

(5.8) 


_ b (~3) + J Ar l/ _ T .. h V-^ as 0, 


iu>o P* - kf 


-^ 3) + , 


ik 


; 2 / 


iu>oP * - kf 


fc(-3) _ 


0, 


XU)Q 1 iojqp* - kf 


• (— 3) 

£>2 = vjJq 


2k 2 , 


(B.9) 

(B. 10) 


( 1 + ¥fp)t[-*> 

\ tuj 0 ) 

2k f 


-(kfn‘P + 1) 


a- 3 ) 

°2 


1 + “ , , ' 
two P ~ kj 

-2P 9( _i ) 6^ 2> -2(l-/i*i , r ( _ 1) )^- 3) = 0, 

(5.11) 


(M*Pr ( _ u - 1)6< 3) - 5 9(-i)t» 


J-2) 


+ (2fc / M*i 5 + l)^~£-^4 3) + iki f pTPb[ 3) 


+ 2 Pkf l -± L b[~ 3) - iki fPbj~ 2) = 0 , 

(5.12) 


+ (2* /M -5 4- i)_^i^-6<- 3) + xfc 2/ M*56'" 3) 

+ 2 Pk f l -^-b[~ 3) - ik 2f Pb { ~ 2) = 0 . 

J lUQ 

(5.13) 

The dispersion relation for ia>o is now obtained 
as follows. From Eqs. (B.6) ~ (B.8), we have 


r(-1) 2fx m P t 1 


,(_3) tfci/ (_ 3) (-3) _ t k 2f (_ 3 ) 

\j\ + 4/j.*P(icJop* + p*Pk 2 { ) j , 

°4 ~ . ‘ "■ t 0 2 j °6 “ • * » °2 i 

IU>0P —kf lU)Qp* — kf 

Jc 2 

l(-3) K f l (-3) 

1 / 

®(-i) = pV P (^o 

+ Pk 2 f ), 

' = — — - 7~rb\ , 

f(-i)(tw 0 p* - fc/) 


(5.2) 

(5.14) 
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Hence, substituting the last of Eqs. (B.14) into the 
difference of Eqs. (B.9) and (B.10), and canceling 
the common factor , we arrive at the dispersion 
relation 


P 

iuop' 


( k f 

- kj \r ( _ i; 


+ 1+^ = 0, (5.15) 


or 


{iw 0 p* -kf)r { _ 1) A’-p*k f r^ 1 )-p*k 2 f =0, (5.16) 

where we note that in going from Eq. (B.15) to 
Eq. (B.16), the extraneous root iujq — kj/p* has 
been introduced into the dispersion relation. It is 
readily shown that for p*P 1, ~ — iujQp*. 

Substituting the latter into Eq. (B.16), we obtain 
a quadratic for Iujq which has the roots kf/A* and 
the extraneous root kf/p *. Thus, for p* P <£ 1, we 
recover the inviscid result iu; o ~ kf/A *, indicating 
instability for A* > 0. 

For the viscous case p*P ~ 0(1), we substi- 
tute the expression for r(_i) given in the first of 
Eqs. (B.2) into Eq. (B.16), and after squaring to 
eliminate the square root and factoring out the ex- 
traneous factor ( iujQp * — kf) noted above, we are left 
with the dispersion relation 

a 2 + Kp'Pkf - 1 )A; - 2p*]a + p*(A; + p') 

-p*Pk f A;(A; + 2p-) = o f 
(5.17) 

where a = iujqp* A* v jkf . Solving for a and choosing 
the nonextraneous root (the one that collapses to 
the inviscid result in the inviscid limit), we finally 
obtain 


Appendix C. 

Effect of Viscosity on the Hydrodynamic 
Pulsating Stability Boundary for k 1 


The existence of the pulsating stability bound- 
ary (51) obtained from the inviscid analysis is highly 
sensitive to viscous effects. Indeed, it is easily shown 
that, in the far outer wavenumber region where vis- 
cous effects are first felt, a weakly nonzero value 
of the liquid-phase Prandtl number 5, while pro- 
ducing no leading-order effect on the inviscid cellu- 
lar boundary, is sufficient to eliminate the pulsating 
type of hydrodynamic instability. Thus, we retain 
the scalings (52) and (62), but consider the weakly 
viscous case 


P = P'yfi, p,P = prp’<? /2 (C.l) 

where the choice of scaling (C.l) permits the same 
scalings for iu; and those coefficients 6* (i = 1, 2, 4, 
6, 8) that apply in the inviscid limit. Thus, based on 
the inviscid analysis of the dispersion relation (49), 
we seek corresponding roots iuj for the weakly vis- 
cous problem in the far outer wavenumber region 
k = kf /e in the form 

iuj ~ e~ 3 / 2 (iuJo + iuj\t 1 ^ 2 H ) , (C.2) 

where, in contrast to the expansion (B.l) in the re- 
gion A* > 0, this form of the expansion will describe 
both the cellular and pulsating hydrodynamic insta- 
bilities in the region A* < 0. 

In the far outer wavenumber region k = kf/e, 
the expansions for r and q are given by 



, kf { A;(i-p'Pk f )-n(A;,k f ), a;> o 
+ 2 ^A;'\A;(i-p*Pkf) + n(A;,k J ), a;< o, 

n(A;,k f ) = J (A;) 2 (i + m *^/) 2 + • 

(5.18) 

It is clear from Eq. (B.16) that for sufficiently small 
p*P and/or |A*|, iuo < 0 for A* < 0. We also 
observe that iwo ~ kf/A * as kf — ► 0, indicating 
that the inviscid result is recovered as the outer 
wavenumber regime is approached. Although this 
root can turn negative for sufficiently large A* (in- 
deed, iuj o = 0 for 2 p*Pk f A* p /p* = 1 - 2p*Pk f + 

yj 1 + 4/i* 2 5 2 A;2 > 0) , this region is of little further 

interest since hydrodynamic instability first sets in 
for negative values of A*. 


r ~ r*(_x/ 2 ) e 1//2 + r o e ° + • * * ; 

r(_i/ 2 ) = p* ~ p*P*k 2 , r 0 = -ivip * , 

(C.3) 

q ~ <?(- 1 ) 6" 1 + <?(_ i/ 2 ) e ~ 1/2 d i 

<?(-!) = {iu 0 /P* +k 2 f ) 1/2 . 

Q(-i/ 2 ) = 1 (iui/P*) (tWo/5* + k 2 y 1/2 , 

(C.4) 

and appropriate expansions for the coefficients 6* are 
written as 


6 i ~ 6|- 2) e - 2 + 6(- 3/2) e - 3 / 2 + - - - , 

(C.5) 

i= 1 , 2 , 8 , 

6 i ~ 6 |- 3 / V 3 / 2 + 6 |- l ) e - 1 + ---, 

(C.6) 

i = 3, 4, 5, 6, 7 . 
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Substituting these expansions into Eqs. (40) - (48) 
and equating coefficients of like powers of we ob- 
tain 

ik i,/ 6^ -3/2) + ik 2 j 4 _3/2) + <7(_ i) b { 7 ~ 3/2) = 0 , 

(C.7) 

ik u b[~ 3/2) +ik 2if b { 6 3/2) +r { -i /2) bi 2) =0, 

(C. 8) 

b[~ 2) = -fc//p* 

(-3/2) 


ifc 2,/ [4 3/2> - 6 4 

-**!,/ [4 


,(-3/2) 


5 -3/2)_ 6 (-3/2)] =| 


(C.9) 

(C.10) 


- ^-b\-‘> = IOJ o , 




-^- 2) + (p* + ^:)6'- 2) = o, 


i + 


2k) P* 

iw 0 


b[~ 2) +bi- 2) 


(C. 11) 
(C.12) 

(CM3) 


-2P* g( _i)4‘ 3/2) -2<>l 


(- 2 ) _ 
'8 ~ 


7 k), 


ifc l./,(-2) _ .,(-3/2) _ ifclj 






P* 

(C.14) 


+ 2P - k _ P- jt! , (,<-*» . % , 

(C.15) 

where Eq. (C.10) was obtained from the leading- 
order difference of Eq. (42) multiplied by ik^f and 
Eq. (43) multiplied by ik\j } and Eq. (C.12) was 
obtained from the leading-order difference of Eqs. 
(44) and (45). 

Equations (C.7) - (C.15) above constitute nine 
equations for the nine leading-order quantities b\ , 
6^~ 2) , b { 3 Z 7 /2 \ b ( a ~ 2) and iu> 0 . Thus, from Eqs. (C.9), 
(C.12) and (C.13), we obtain 


( 1+ m :w* f _ 1 p W 

\ IU)Q J 1 p* 


(-3/2) 

7 




(C.16) 

while the sum of Eq. (C.14) multiplied by ik\j 
and Eq. (C.15) multiplied by along with Eqs. 
(C.7) and (C.9), gives 

[??-„ + fc ?R J/2> -f£‘'.- 2> = o. < C17 > 


Hence, using the definition of Eqs. (C.ll) and 
(C.17) determine b[~ 2 ^ and b?~ 3 ^ as 

b (_ 2) = — -j— — (2P*fc? + iu) o) , 6y _3/2) = -2 P*k) . 

kj J 

( c . 18) 

Substituting the results (C.17) and (C.18) into Eq. 
(C.16), we finally obtain the leading-order dispersion 
relation for iu o as 


(iu 0 ) 2 4- 4 P' k } iu)Q - 4P*kj fa oP* 4- k 2 P* 2 = 


2k 2 


- A;) - 4 p-V, , 

(C. 19) 

where A* is given by the expression (55) for the 
leading-order inviscid cellular boundary A in the 
far outer wavenumber region, and the square root, 
which arises from the definition in (C.4) of 
denotes the principal square root. 

In the inviscid limit P* = 0, we recover the 
leading-order results for the inviscid case. In par- 
ticular, Eq. (C.19) reduces to {iuo) 2 = 2A^(A* - 
Ap)/p* 2 i and hence iu 0 > 0 (unstable) for A p * > A* 
and is pure imaginary for A* < A*. Stability in the 
latter region is thus determined by the coefficient 
iu\ of the next term in the expansion (C.2), which, 
from either a continuation of the present calculation 
or from the direct expansion of the inviscid disper- 
sion relation (49), is given by iu i = (A* 2 ~p* 2 )/p* 2 . 
Consequently, the region —p* < A* < A* is stable, 
while the region A* < — p* is unstable, and identical 
results are obtained in the inner and outer wavenum- 
ber regimes. Hence, we recover both the steady (cel- 
lular) stability boundary A* = A* (on which iu is 
identically zero) and the nonsteady (pulsating) sta- 
bility boundary A* = - p * (on which only the real 
part of iu is zero) shown in Figures 2 and 3. 

For P* > 0, it is readily seen from Eq. (C.19) 
that the root iu o = 0 still corresponds to the invis- 
cid cellular boundary A* = A*, so that the cellular 
boundary is unchanged from its inviscid limit in the 
weakly viscous regime described by (C.l). On the 
other hand, it can be shown that the entire region 


a; < 


Ap is now stable. 


This is most easily seen 


by considering the regime 0 < P* 1 and seeking 
solutions iu 0 of Eq. (C.19) in the expanded form 


lUQ 


' iao 4- ia\P* 4- ic*2P* 2 4- 


(C.20) 


Then, at leading order with respect to the above ex- 
pansion, we obtain the previous inviscid result (cor- 
responding to P* =0), 


2k 2 

(*«o) 2 “ RW - A ;) , 


(C.21) 
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which implies instability only for A* > A*, and neu- 
tral stability (Re{iao} = 0) otherwise. At the next 
order, however, we obtain 

zc*o(iai 4- 2 k 2 f ) = 0 , (C. 22) 

which implies stability (ia i = -2 kj < 0) in the re- 
gion A* < A*. Thus, in the A p regime considered 
here, even very weak (P* < 1) viscosity is suffi- 
cient to destroy the pulsating type of instability, at 
least in the far outer wavenumber region where vis- 
cous effects are first felt, whereas the cellular stabil- 
ity boundary is clearly unaffected even when P* ~ 
0(1). Indeed, the latter is modified only for 0(1) 
values of the unsealed Prandtl number P, as de- 
scribed in the main body of the text. 


Appendix D. 

The Pulsating Stability Boundary for 
Small Liquid Prandtl Numbers 


For small liquid viscosities such that P = Pe, 
and the same scalings (52) and (62) as in the P ~ 
0(1) case analyzed in the main body of the text, 
the effect of viscosity on the cellular boundary dis- 
appears at leading order [P — ► 0 in Eqs. (100)], while 
the effect of viscosity on the pulsating boundary is 
only significant in the far outer wavenumber regime. 
In that case, the appropriate expansions analogous 
to Eqs. (101) - (106) are given by 

iu ~ e 3/2 (iujo 4- ioJie 1/4 + itJoe 1/2 4 ) , (Z?.l) 


r ~ r^i/ 2 )C 1 ^ 2 + 0(e 1 ^ 4 ) , ^(— 1 / 2 ) = — > 

(D. 2) 

Q ~ 9(-5/4)C -5 / 4 + 0(e -1 ) , <J(-5/4) = (iwo/P) / , 

(D. 3) 

b t = b[~ 2) t- 2 + &<- 7/4 >e" 7 / 4 + b^t-V 2 + ■ • ■ , 

t = 1 , 2, 8 , 

{DA) 

bi = 6 (-3/2) e -3/2 +6 (-5/4) £ -S/4 + . , = 3> 4 , 5, 6 , 

(D. 5) 

bi = 6<- 5/4) €- 5 / 4 +6[- 1) £- 1 +• • • , 1 = 7. (D. 6) 


Substituting these expansions into Eqs. (40) - (48) 
and equating terms corresponding to like powers of 
e then gives, as previously, a sequence of equations 
for the recursive determination of the coefficients in 
the above expansions. Similar to the calculation for 


P ~ 0(1) in the outer wavenumber regime (Section 
6), we obtain in this case that 

W 2 = 4K'4 iwi = 0 , 

9 r 2 „ (D.7) 

iuj 2 = kj I (^r) - 2Pk/ - 1 , 

where A* = {p* /2)(p*1fkf - 1) is the inviscid cel- 
lular boundary in the far outer wavenumber regime 
given by Eq. (55). The first of Eqs. (D.7) thus re- 
covers the cellular stability boundary, but since iwo 
is purely imaginary for A* < A*, stability in that 
region is determined by the real part of the next 
nontrivial coefficient in the expansion (D.l). Thus, 
setting iu>2 = 0 in the last of Eqs. (D.7), the pulsat- 
ing stability boundary in the far outer wavenumber 
regime is given by 

A; = -p*y/l + 2Pk f , (D. 8) 

which, in the limit k/ — * 0, matches with the leading- 
order pulsating boundary A* = —p* in the outer 
wavenumber region, which is unaffected by viscosity 
to this order of approximation. Thus, Eq. (D.8), 
which in terms of k is given by A* = —p*( 1 + 

2 Pek) ' , is valid for arbitrary wavenumbers. Writ- 
ing P in terms of its unsealed counterpart P, this 
expression becomes identical to Eq. (129), which 
thus remains valid in the limit of small P. 
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